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Abstract 

The knowledge of stellar evolution is evolving quickly thanks to an increased 
number of opportunities to scrutinize the stellar internal plasma properties by 
stellar seismology and by ID and 3D simulations. These new tools help us to 
introduce the internal dynamical phenomena in stellar modeling. A proper 
inclusion of these processes supposes a real confidence in the microscopic 
physics used, partly checked by solar or stellar acoustic modes. In the present 
paper we first recall which fundamental physics has been recently verified 
by helioseismology. Then we recall that opacity is an important ingredient 
of the secular evolution of stars and we point out why it is necessary to 
measure absorption coefficients and degrees of ionization in the laboratory 
for some well identified astrophysical conditions. We examine two specific 
experimental conditions which are accessible to large laser facilities and are 
suitable to solve some interesting questions of the stellar community: are 
the solar internal radiative interactions properly estimated and what is the 
proper role of the opacity in the excitation of the non radial modes in the 
envelop of the /3 Cepheids and the Be stars ? At the end of the paper we point 
out the difficulties of the experimental approach that we need to overcome. 

Key words: stellar plasma, opacity coefficients, reaction rate, absorption 
spectra 



1. The context 

The stellar structure is currently determined by the resolution of four 
equations which govern the hydrostatic equilibrium, the conservation of mass, 
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the energetic equilibrium and the energy transport dominated by radiation 
or convection. In order to describe the secular temporal evolution of stars, 
one needs to add to these coupled equations the temporal evolution of the 
composition through nuclear reaction rates and diffusion or turbulent effects 
PQ. In fact the resolution of these equations supposes a good description of 
plasma properties on a very broad range of temperature (typically several 
thousand to billion degrees) and density (10~ 7 to 1000 g/cm 3 ) in term of nu- 
clear interaction, equation of state and photon interaction with the different 
elements of the stellar plasma. Today helio and astero seismology allow to 
check the hypotheses of such a calculation and new indicators are extracted 
from helioseismology to progress towards a more dynamical view of stars 
(2} IB]. Such approach is now extended to other stars with a better under- 
standing of massive stars (U |5]. This crucial effort gives us the potential to 
solve new problems in our knowledge of stars: 

- the development of a unified theory of stellar evolution will naturally link 
short time stages dominated by the dynamics (stellar formation, young stars 
modeling or supernovae explosion) and long stages dominated by slow mi- 
croscopic evolution but where dynamical effects continue to be present, 

- the complete MHD description of the Sun will connect internal and external 
dynamics. This will help to establish the origin (s) of the solar activity and its 
prediction for the next century with its consequent impact on the terrestrial 
environment. Such approach will also be useful for understanding stellar 
winds in massive stars which turn quickly, are deformed and are magnetically 
active. 

Consequently, we are currently in a period of transition where we evolve 
from a qualitative description of stars to an entirely quantitative one. On one 
hand, one can say that some plasma properties are definitively established 
(section 2). On the other hand, internal rotation and magnetic field induce 
transport of momentum and chemicals which are still difficult to estimate 
properly due to our poor knowledge of the initial conditions and of the in- 
creased number of transport coefficients in the equations. So there is a real 
benefit to isolate some phenomena and go to the laboratory for a verification 
or an improvement of the theoretical predictions. 

Three decades ago, stellar evolution drove measurements of nuclear cross 
sections in the laboratory. This activity is now oriented toward the ultimate 
stages of evolution. Today large laser facilities have become very crucial tools 
to check the absorption coefficients in some specific conditions, this informa- 
tion is fundamental to verifying the transport of energy inside the stars by 



2 



photon interaction and to improve our description of the stellar envelops. 
Section 3 will examine how the absorption coefficients are introduced in the 
stellar equations, we deduce from them two examples where opacity coef- 
ficients must be controlled. The astrophysical need leads to a campaign of 
absorption coefficient measurements as fundamental checks of atomic physics 
in stellar plasma conditions like it has been done for nuclear cross sections 
in the last two decades. 

Section 4 will focus on the specific needs and on the natural limitations of 
the experimental side that we must reduce if we want that these experiments 
will contribute, among others, to the renewal of stellar evolution. 

2. Stellar plasma properties checked by helioseismology 

Helio(astero) seismology is a discipline exploring the standing waves which 
are permanently excited by the superficial or deep convection of the Sun or 
stars. It has revealed interesting internal properties of the solar plasma. 

Acoustic waves propagate inside the whole Sun and their travel time 
depends on the local sound speed. They generate very slightly motions in 
the atmosphere of the Sun which are detectable. We can formely treat this 
information through a perturbative theory because these perturbations are 
small and the Sun is reasonably spherical. 

The Sun (or a solar like star) is a self-graviting sphere of compressible 
gas; it oscillates around its equilibrium state with a period of about 5-min. 
These oscillations are interpreted as a superposition of waves propagating 
inside the star (acting as a resonant cavity), and forming standing waves: the 
eigenmodes of vibration. By projecting these modes onto spherical harmonics 
Y/ 11 , we write any scalar perturbations as (in the case of the pressure p ) [6]: 



£{r,6,ip,t) = lfa{r),fa{r)-.,fa{r)^-— y™(9, cp) expiu^t (1) 



where fa = l/(u 2 r)[p /p + $'] is the horizontal displacement, $ the gravita- 
tional potential perturbation, u n j )m the eigenfrequency and p the gas density. 
The quantum numbers n, I, m are respectively the radial order (number of 
nodes along the radius), the degree (the total horizontal wave number at the 



p'(r, 9, ip, t) = p'(r)Y i m (6', ip) exp iu n ^ m t 



and the vector displacement £ as 
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Figure 1: Eigenfunctions of an acoustic mode (i =0 radial mode, n=23), and of a gravity 
mode [l =2, n=10). These modes are now detected in different stellar sites. 



surface is k^ ~ L/R Q , with L = \Jl{l + 1)) and the azimuthal order (number 
of nodes along the equator with \m\ < I). Restricting the phenomenon to 
adiabatic oscillations within the Cowling approximation (<3> neglected) and 
considering only small radial wavelengths compared to R®, we reduce a 4th- 
order system equations to a second-order wave equations, with the following 
dispersion relation [5] 



K- c2 




2 2 
+ U n,Lm ~ U c 



(2) 



where the squared length of the wave vector is written as the sum of a radial 
and a horizontal component: \k\ = k^ + k 2 ., k\ = Ff /c 2 s is the horizontal wave 
number, Ff = L 2 c 2 s /r 2 the Lamb frequency, iV 2 = g[l/Tid]np/dr — dlap/dr] 
the Brunt- Vaisala frequency, uj 2 c = c 2 (l — 2dH p /dr)/AH 2 the acoustic cut- 
off frequency (~ 5.8 mHz), H^ 1 = —d\np/dr the density scale height, Ti 
the adiabatic exponent and c 2 = Tip/p the sound speed. The oscillatory 
solutions of the wave equation define two types of waves, i.e acoustic ones 
(with uj n ^m > N,Fe) and gravity ones (with uj n ^ m < N,F(). 

Figure 1 shows that the acoustic modes have their maximum amplitude 
at the surface (left side) but the low £ modes penetrate down to the core. 
The gravity modes are excellent probes of the solar (stellar) core and are 
evanescent at the surface with rather small amplitude (right side). A lot 
of solar acoustic modes and some gravity modes have been detected. So a 
refined analysis of their properties leads to a stratified analysis of the solar 
internal structure from the surface down to the inner 6% of the radius thanks 
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to a very high accuracy (better than 10~ 4 ) extraction of the sound speed [7]. 
Figure 35 of pQ shows the sensitivity of the sound speed profile through 
the ratio T/fi where \i is the mean molecular weight to physical processes. 
The inflexion of the profile due to the partial ionisation of helium or due 
to the change of energy transport between radiation and convection (see 
Figure 46 of pQ) has led to the determination of the solar helium content 
and of the position of the basis of the convective zone. Recently, two other 
important properties of the dense central plasma have been checked from 
the sound speed profile of the nuclear core: the maxwellian distribution of 
the reactant particle velocities and the absolute value of the proton proton 
reaction rate [8] in order to properly predict the solar neutrinos and also the 
central temperature. 

The sensitivity of the sound speed to the different ingredients of the solar 
model has been studied in [9] . This observable is particularly sensitive to the 
detailed solar composition through the opacity coefficients and one shows 
that a l%variation of the opacity coefficient produces an 0.1% variation of 
the sound speed which is known with an accuracy of some 10~ 5 [7J. The 
use of the updated determination of the CNO superficial composition has 
led to a large discrepancy between the observed sound speed and the sound 
speed extracted from the standard solar model. In fact this discrepancy is 
about 15 to 26 times the vertical uncertainty as shown on Figure 2 [TU]. The 
two recent estimates of the solar composition are now confirmed by a third 
estimate and the discrepancy converges towards the curve corresponding to 
what we call here the Holweger composition. 

The interpretation of this discrepancy is not yet understood. A lot of 
potential interpretations have been studied, see for example [TT] and the 
review of [12] which gives a summary of the present hypotheses. The origin 
could be an insufficient knowledge of the composition inside the radiative 
zone or of the related opacity coefficients. The effect is about 20-30% or 10% 
on the mean Rosseland values depending if the problem is connected to only 
one element or to most of the heavy elements (Z > 2). Another explanation 
could be that the solar model is incomplete because the dynamical effects 
are not yet introduced and could provide some transport of the different 
species. Helioseismology cannot solve this problem due to the connection 
between abundance and opacity within the solar model. Consequently this 
important discrepancy encourages us to study some experimental conditions 
in the laboratory and to measure appropriate absorption spectra to validate 
or improve the complex theoretical estimates. 



5 



0.03 r 



0.02 



0.01 7 




0.0 0.2 0.4 0.6 0. 

r/R G 



8 1.0 



Q, 




1.0 



Figure 2: Squared sound speed and density difference deduced from the SoHO observations 
and the solar model using 2 updated compositions: Holweger (- . - .-) and Anders (full 
line + error bars) or with the seismic model. See details in [TO]. The recent reanalysis of 
the solar oxygen and nitrogen estimate favors a difference similar to the one called here 
Holweger. 
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3. The stellar opacities 



Stellar opacities partly define the longevity of stars. In fact the lifetime 
of the stars depends on the stellar mass and the heavy element contribu- 
tion to the total opacity coefficient while the stellar composition is generally 
dominated by hydrogen and helium. In fact, except below the surface, the 
light elements are generally totally ionized so their contribution to the ab- 
sorption coefficients is easier to estimate. Inside stars, the energy transfer in 
the radiative zones follows the equation: 

^ = zi^im where jl = r^ -^AT)]^^^ 

dr Aac T 3 4vrr 2 k k f °° V ; 

kr is the Rosseland opacity at the position r and for a mixture of species 
Xj. The weighting factor gives the greatest weight to the region around 
v = AkT/h. L, T, p are the local luminosity, temperature and density of the 
plasma at the position r. In these conditions, a species which is not totally 
ionized and has a large Z (number of protons) value may play an important 
role [13J. K ua is the term of absorption, it includes free-free, bound-free and 
bound-bound processes and k vs is the scattering term. 

3.1. Sun and solar-like stars 

The calculation of the opacity coefficients is complex for solar like stars 
because the plasma is generally partially ionized so the astrophysical com- 
munity has greatly benefitted from the work of the Los Alamos team (LAOL 
tables) [16j and more recently from the OPAL team [151 E] who has deliv- 
ered tables for a large range of temperature, density and composition since 
1996. A third academic team [20] has first worked on the opacities of the 
stellar envelop and extended their work to all stellar conditions in providing 
the OP tables. This team provides also the corresponding spectra [2T] . 

Figure 3 shows the respective role of different elements in the internal solar 
conditions for the two used solar composition and OPAL opacities [TU t [18], H~9] . 
the same estimate using the OP opacities leads to slight differences mainly 
for the oxygen contribution which varies by a little more than 10% with an 
impact of the total Rosseland coefficient of about 4% near the convective 
basis [22]. The large contribution of iron, neon and oxygen in the radiative 
zone (below 0.71 R ) is due to their bound-bound contributions. Silicon also 
plays a role at a level of 10% below 10 7 K. It is interesting to note that iron 
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Figure 3: Relative contribution to the total opacity of the most abundant species which 
play an important role in the opacity of the radiative zone after hydrogen and helium. This 
estimate uses the OPAL opacity calculations for the Grevesse and Asplund composition 
differing by the content of C, N, O. 

has a fraction number of 5 10 -5 and is always partly ionized even in the center 
of the Sun, its contribution to the total opacity is at a level of 20% in most of 
the radiative zone. Oxygen becomes partly ionized at 0.4 R Q and plays the 
major role to limit the radiative transfer, around 0.71 R Q . The increase of its 
opacity contribution triggers the convection which is the dominant process 
in the transport of energy in the outer layers. 

For higher stellar masses the central temperature increases and all the 
species are totally ionized (see for example figure 2 of [T3] or 31 of [TJ). In 
these stars the core is convective due to the amount of energy produced by 
the role of CNO, so the lifetime of these stars is mainly driven by the nuclear 
reaction rates rather than by opacity coefficients. Figure 4 illustrates this 
point in showing the role of iron (for a solar composition) at the limit of the 
convective core for different stellar masses. 

3.2. The microscopic diffusion of heavy elements and the k mechanism 

We now examine two other important roles played by the opacity coeffi- 
cients. The life of solar-like stars may reach several billions of years. During 
this long period of time, the different species migrate slightly towards the 
stellar core due to the gravitational settling of the elements heavier than hy- 
drogen. This effect is extremely slow, resulting for the solar case in 10-15% 
depletion at the stellar surface after billion years well measured for helium. 
In stars a little more massive stars than the Sun, this effect is partly compen- 
sated by the levitation of some species due to radiative pressure [231 El] • To 
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Figure 4: Relative role of iron, oxygen and silicium (for a solar composition) at the limit 
of the convective core for different stellar masses. 



estimate the respective role of gravity and pressure, one needs to calculate 
the radiative acceleration for each species k: 

9rad{k) = T^Wf Rl{k) Wh6re 7(k) = / ^al) IAV (4) 

where Fr is the radiative flux. In such a calculation the detailed of the 
opacity spectrum is needed and f u is a function of the frequency and of 
the temperature. It seems that there is still some large discrepancies in 
the radiative acceleration between the two opacity tables available for the 
international stellar community, the OPAL and the OP tables [21] • From 
carbon to iron one may notice some differences up to 50% between 1.2 eV 
and 100 eV, in particular for iron [22J, the origin of the difference is not yet 
understood and could be investigated by some laboratory experiment. 

Another important process, called the k mechanism, is at the origin of 
the pulsation of intermediate stars. There is a renewal of interest for that 
mechanism since the launch of COROT (2006) which is the first satellite ded- 
icated to the development of asteroseismology [26] and search for exoplanets. 
This mechanism is responsible for the pulsational instability of stars between 
1.5 to 20 M and justifies some explanation. 
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Figure 5: Relative role to the total opacity of the different elements for 1.85 and 12 M Q 
through the whole star. The profile of the total Rosseland opacity of the mixture (in AU) 
is shown by a thin black solid line. The important role of iron is clearly visible in the 
external radiative zone and creates the k mechanism when these layers extend or contract 
with time [35] • (See the color version of the figure in the electronic file). 



In the external optically thick layers, the diffusion approximation is jus- 
tified and: 

1 d5L r 8L T dr dn A .5T St. 

The first term describes the radiative dissipation and in fact stabilizes the 
star, the term ^ describes the direct influence of the temperature variation 
on the luminosity, this effect is generally small, the last term contributes 
to the instability because in case of the compression the radiating area is 
reduced. The k mechanism will work if the opacity varies more quickly 
with radius than the other terms. Two kinds of excitation are possible: the 
opacity bump is due to the partial ionization of He (He + to He ++ ) at log T 
approximatively 4.5 and excites the evolved cepheids, RR Lyrae (low mass 
evolved stars) and the 5 Scuti stars (1.5-2.5 M© main sequence stars) or 
the opacity bump is due to the M shell of Fe at log T approximatively 5.2 
and excites the (3 Cepheids stars (7-20 M© main sequence stars), SPB (3-9 
M© main sequence stars) and sdB stars (helium core, 0.5-1.4 M©). In a /3 
Cepheid star of 8 MQ the total opacity reaches the value of 6.95 cm 2 g _1 due 
to the bound-bound process associated to the partial ionization of iron at 
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T= 195800 K (leV= 11600 K) and for a density of 6.6 10~ 7 g/ cm3 - We show 
on Figure 5 the contribution of iron and some other elements in the external 
radiative zone of different stars from 1.8 to 12 M Q [25] . 

The Cepheid stars are known to be used to estimate the stellar distances. 
They have been better understood after a revision by a factor 2 of the opacity 
for species of Z> 2 confirmed by the first measurements of absorption coef- 
ficients. The (3 Cepheids stars are now of first interest due to the beautiful 
new asteroseismic observations and the correlation between the frequencies 
of the excited modes and the mass of these stars for a given age. 

These stars are now used to estimate the age of several young open clus- 
ters [27J. Nevertheless the labeling of the detected frequencies is not an easy 
task so complementary observations or measurements are needed. The dif- 
ference between OP and OPAL tables leads to some interesting differences 
in the stellar surface temperature and it is presently believed in the seismic 
community that the use of OP opacities seems to help the interpretation of 
the observed frequencies [2"g] . 

4. From theoretical estimates to laboratory experiments 

We deduce from the above stellar situation two kinds of conditions which 
will benefit from the present and coming large laser facilities. 

The interest to describe better the physics of the stellar envelope is jus- 
tified by the study of more and more pulsating stars (COROT, KEPLER) 
with the objectives to put constraints on the corresponding stellar interiors. 
This fact leads to study relatively low temperature (about 20 eV) and very 
low density. Of course these densities are too low to be considered presently 
and one needs to search equivalent situation at higher density. This kind 
of investigation could be done on medium size lasers like the LULI2000 in 
France and also GEKKO XII and VULCAN, an interesting case is shown in 
[29] when we compare different neighbors of relatively high Z at 20 eV and 3 
mg/cm 3 . For the astrophysical purpose we will enlarge the spectrum domain 
to XUV between 10 to 100 eV to cover the most important stellar region 
where the Rosseland mean value is maximal. For these conditions, OP and 
OPAL spectra present some different behaviors and we are preparing the best 
experimental conditions to resolve this discrepancy and improve the stellar 
interpretation of the seismic data [30] . 

The second case is connected to opacity coefficients for the internal struc- 
ture of radiative zones of solar-like stars. The temperature is generally greater 
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than 100 eV and the density greater than fraction of solid density. If we want 
to investigate different cases, one needs to work on coming larger facilities 
(LIL +PETAL, OMEGA EP, FIREX II, LMJ, NIF) because the difficulty 
is mainly to reach the equivalent high density at LTE with a precise diag- 
nostic. In these cases, on needs to use an X-ray spectrometer to measure 
the frequency band necessary for the astrophysical purpose or to find some 
equivalence to check the related plasma properties. Let consider this second 
case in more details. 

4-1. Degree of ionization of the different species in the solar plasma 

The composition of the solar plasma is by definition different than the 
plasma produced in laboratory. So we estimate here the differences. Let us 
first reduce the solar plasma to a 4 components plasma formed by H, He, O 
and Fe, with the respective fraction numbers are He = 0.08, = 5 10 -4 , Fe = 
2.8 10 5 , hydrogen accounting for the remainder (following approximatively 
the solar composition). Clearly in the solar plasma, the free electrons come 
from hydrogen so it is interesting to calculate the degree of ionization of Fe 
and O in considering pure iron or oxygen or in considering the simplified 
solar composition. The calculations have been made for two conditions: at 
the basis of the convective zone where the temperature is about 193 eV and 
the density 0.2 g/cm 3 and for the iron at the maximum of the iron opacity 
peak around 400 eV and a density of 2.5 g/cm 3 . 

Table 1 shows the averaged degree of ionization for iron and oxygen ob- 
tained with different codes: the FLYCHK code which provides ionization and 
population distributions of plasmas using a schematic atomic structures and 
scaled hydrogenic cross- sections [3T] that we use for one or two components 
(H, Fe or H, O), a mixture model using Thomas Fermi averaged ionization 
[32] which allows us to consider the 4 components and also the OP code 
where we consider only pure Fe and O. One notices that the choice of the 
mixture does not significantly change the degree of ionization when T and 
N e (electron number) are fixed but that this degree of ionization depends on 
the hypotheses of the model used (mainly for iron). The density of the cor- 
responding pure iron measurement at 193 eV for the same averaged degree 
of ionization can be obtained in replacing the density of the plasma by p#p 
where £* is the ionization degree. In the discussed case the density of the 
pure iron must be = 0.6 g/cm 3 instead of the value of 0.2 g/cm 3 of the solar 
plasma. Of course an opacity experiment of a mixture of elements could be 
also another interesting experiment to perform. 
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Tabic 1: Average degree of ionisation for iron and oxygen for different density and 
temperature and mixture conditions obtained by three calculations: FLYCHK (F), OP 
(OP) and a mixture model using Thomas Fermi (TF). We consider different cases: pure 
iron or oxygen, the same element in a mixture of H+ Fe (or O) and the same element in 
a mixture of 4 elements. 
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Figure 6: Comparison of the individual iron and oxygen spectra for different conditions 
of temperature and density around the BCZ solar conditions in using the OP tables. (See 
the color version of the figure in the electronic version). 



4-2. Sensitivity of the absorption spectra to the temperature and density con- 
ditions 

The Z pinch facility of Sandia National Laboratory has been used for 
the first comparison between experiment and calculations [33]. The exper- 
iment has been done at T e = 156 ± 6 eV and N e = 6.9 ± 1.7 10 21 cm -3 . 
These thermodynamic conditions seem not so far from the solar conditions 
discussed in table 2 even it is always difficult to determine the real conditions 
of measurements [34] . 

The experimental difficulties in getting a detailed comparison with the 
spectra used in astrophysics are numerous. First one needs to determine 
without ambiguity the proper temperature and density of the experiment, 
then one needs to obtain the proper level of transmission to look to the 
absolute value of the detailed spectra. We illustrate on Figure 5 how the 
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iron and oxygen spectra evolve for the OP calculations with temperature 
and density in a range probably larger than the uncertainty we may have for 
these two values. Doing so we illustrate directly the complexity of the iron 
spectra in comparison with the oxygen spectra. Of course the interesting 
spectrum range is around 800 eV which corresponds to the maximum of the 
mean Rosseland value. Although dominant peaks present small variations, 
peak to valley ratio is clearly different. They are sensitive to temperature 
variation of 10% and yields to a variation of 6% in Rosseland opacity in the 
case of iron, the spectra show clearly the difference in ionization numbers. 
The effect on the Rosseland mean value for a change of density by a factor 3 
is between 8 and 13%. 

So several measurements of absorption spectra obtained on different fa- 
cilities will be useful. They require a good determination of the plasma tem- 
perature and density which can be delivered by petawatt lasers and a well 
calibrated transmission. These experiments will help to disentangle the ori- 
gins of the observed differences between the calculations from OP and OPAL 
for some astrophysical conditions. They will also give strong confidence in 
these very complex atomic calculations and will contribute to build a more 
sophisticated view of stars which must include microscopic and macroscopic 
phenomena to reproduce all the seismic observables: sound speed, superficial 
helium content, internal rotation, circulations, magnetic field. 
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